Chapter 1 



Time-Frequency Concepts 



Time-frequency signal analysis and processing (TFSAP) concerns the analysis and 
processing of signals with time-varying frequency content. Such signals are best 
represented by a time- frequency distribution (TFD), which is intended to show 
how the energy of the signal is distributed over the two-dimensional time-frequency 
space. Processing of the signal may then exploit the features produced by the 
concentration of signal energy in two dimensions (time and frequency) instead of 
only one (time or frequency). 

The first chapter begins the introductory tutorial which constitutes Part I of 
the book. This tutorial updates the one given in [1] by including recent advances, 
refining terminology, and simplifying both the presentations of concepts and for- 
mulations of methodologies. Reading the three chapters of Part I will facilitate the 
understanding of the later chapters. 

The three sections of Chapter 1 present the key concepts needed to formu- 
late time-frequency methods. The first (Section 1.1) explains why time-frequency 
methods are preferred for a wide range of applications in which the signals have 
time- varying characteristics or multiple components. Section 1.2 provides the signal 
models and formulations needed to describe temporal and spectral characteristics 
of signals in the time-frequency domain. It defines such basic concepts as analytic 
signals, the Hilbert transform, the bandwidth-duration product and asymptotic 
signals. Section 1.3 defines the key quantities related to time- frequency methods, 
including the instantaneous frequency (IF), time-delay (TD) and group delay (GD). 



0 Author: Boualem Boashash, Signal Processing Research Centre, Queensland University 
of Technology, Brisbane, Australia. Reviewers: K. Abed-Meraim, A. Beghdadi, M. Mesbah, 
G. Putland and V. Sucic. 
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Chapter 1: Time-Frequency Concepts 



1.1 The Need for a Time- Frequency Distribution (TFD) 

The two classical representations of a signal are the time-domain representation s(t) 
and the frequency-domain representation S(f). In both forms, the variables t and / 
are treated as mutually exclusive: to obtain a representation in terms of one variable, 
the other variable is “integrated out”. Consequently, each classical representation 
of the signal is non-localized w.r.t. the excluded variable; that is, the frequency 
representation is essentially averaged over the values of the time representation at 
all times, and the time representation is essentially averaged over the values of the 
frequency representation at all frequencies. 

In the time- frequency distribution, denoted by the variables t and / 

are not mutually exclusive, but are present together. The TFD representation is 
localized in t and f. 

1.1.1 Representation of Three Real-Life Signals 

The usefulness of representing a signal as a function of both time and frequency is 
illustrated by considering three signals of practical importance: 

1. Sinusoidal FM signal: Monophonic television sound, like monophonic FM 
radio, is transmitted on a frequency-modulated carrier. If the audio signal is 
a pure tone of frequency / m (the modulating frequency), then the frequency 
of the carrier is of the form 

fi(t) = f c + /dC0s[27r/ m t + (f>] (1.1.1) 

where t is time, fi(t) is the frequency modulation law (FM law), f c is the 
mean (or “center”) carrier frequency, /< j is the peak frequency deviation, and 
< f> allows for the phase of the modulating signal. The amplitude of the carrier 
is constant. 

2. Linear FM signal: Consider a sinusoidal signal of total duration T, with 
constant amplitude, whose frequency increases from fo to /o + 13 at a constant 
rate a = B/T. If the origin of time is chosen so that the signal begins at t = 0, 
the FM law may be written 

fi{t) = fo + at ; 0 <t<T. (1.1.2) 

In an electronics laboratory, such a signal is called a linear frequency sweep, 
and might be used in an automated experiment to measure the frequency 
response of an amplifier or filter. In mineral exploration, a linear FM signal 
is called a chirp or Vibroseis signal, and is used as an acoustic “ping” for 
detecting underground formations [2,3]. 

3. Musical performance: A musical note consists of a number of “components” 
of different frequencies, of which the lowest frequency is called the fundamental 
and the remainder are called overtones [4, p. 270]. These components are 
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Fig. 1.1.1: An example of musical notation [5]. Roughly speaking, the horizontal dimension is time 
and the vertical dimension is frequency. 



present during a specific time interval and may vary in amplitude during that 
interval. In modern musical notation, each note is represented by a “head”. 
The vertical position of the head (together with other information such as 
the clef and key signature) indicates the pitch, i.e. the frequency of the most 
prominent component (usually the fundamental). The horizontal position 
of the head in relation to other symbols indicates the starting time, and the 
duration is specified by the shading of the head, attached bars, dots, stem and 
flags, and tempo markings such as Allegro. The choice of instrument — each 
instrument being characterized by its overtones and their relationships with 
the fundamental — is indicated by using a separate stave for each instrument or 
group of instruments, or a pair of staves for a keyboard instrument. Variations 
in amplitude are indicated by dynamic markings such as mp and crescendo. 
Fig. 1.1.1 illustrates the system. By scanning a set of staves vertically, one can 
see which fundamentals are present on which instruments at any given time. 
By scanning the staves horizontally, one can see the times at which a given 
fundamental is present on a given instrument. 

Each of the three above signals has a time-varying frequency or time-varying 
“frequency content”. Such signals are referred to as non-stationary. 

The three examples described above are comprehensible partly because our sense 
of hearing readily interprets sounds in terms of variations of frequency or “frequency 
content” with time. However, conventional representations of a signal in the time 
domain or frequency domain do not facilitate such interpretation, as shown below. 



1.1.2 Time-Domain Representation 



Any signal can be described naturally as a function of time, which we may write 
s(t). This representation leads immediately to the instantaneous power, given 
by |s(i)| , which shows how the energy of the signal is distributed over time; the 
total signal energy is 



E = 



|s(f)| 2 dt. 



( 1 . 1 . 3 ) 
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But the time-domain description has limitations, as may be seen by applying it to 
the above three examples: 

1. The sinusoidal FM signal whose frequency satisfies Eq. (1.1.1) may be writ- 
ten 

si ( t ) = A cos (^2tt f c t + j ^ sin[27r f m t+<j>] + V’) (1-1.4) 

where A is the amplitude and ip a phase offset; the fraction /d// m is called the 
modulation index and is equal to the peak phase deviation (in radians) from 
2iTf c t. This equation, by itself, does not clearly show how the frequency varies 
with time. If we imagine a graph of s\(t) vs. t , it would give the impression of 
an oscillating frequency, but determining the frequency as a function of time 
from the graph would be a tedious and imprecise exercise. 

2. The linear FM signal whose frequency satisfies Eq. (1.1.2) may be written 



S 2 (t) = A rect 



t-T / 2 
T 



cos (27r[/ 0 t + ft 2 ] + Ip) 



(1.1.5) 



where, again, A is the amplitude and ip a phase offset. The rect function is 
a rectangular pulse of unit height and unit duration, centered on the origin of 
time; that is, 



rect t = 



1 if M < 1/2 ; 

0 otherwise. 



( 1 . 1 . 6 ) 



Hence the rect[. . .] factor in Eq. (1.1.5) is equal to unity for 0 < t < T, and 
zero elsewhere. But again it is not immediately apparent why Eq. (1.1.5) has 
the required FM law. 

The graph of S 2 (t) vs. t is shown on the left side of Fig. 1.1.2(a), for ip = 0, 
T = 64s, fo = 0.1 Hz and a = (3/640) Hzs -1 . Although the graph gives a 
clear impression of a steadily increasing frequency, the exact FM law is not 
clearly and readily visible. 



3. A musical performance can be represented as (for example) an air pressure 
curve at a particular point in space. Each such curve is a time-varying pres- 
sure, and may be converted by a microphone and amplifier into an electrical 
signal of the form S'i(t). Indeed, music is routinely recorded and broadcast 
in this way. However, the function ss(t) is nothing like the form in which 
a composer would write music, or the form in which most musicians would 
prefer to read music for the purpose of performance. Neither is it of much use 
to a recording engineer who wants to remove noise and distortion from an old 
“vintage” recording. Musical waveforms are so complex that a graph of s^pt) 
vs. t would be almost useless to musicians and engineers alike. 



These three examples show that the time-domain representation tends to obscure 
information about frequency, because it assumes that the two variables t and / are 
mutually exclusive. 
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Fig. 1.1.2: The importance of phase: (a) time-domain representation (left) and magnitude spectrum 
(right) of a linear FM signal [Eq. (1.1.5)] with duration 64 seconds, starting frequency 0.1 Hz and 
finishing frequency 0.4 Hz; (b) time-domain representation and magnitude spectrum of another signal 
with the same magnitude spectrum as that in part (a). All plots use a sampling rate of 8 Hz. 



1.1.3 Frequency-Domain Representation 

Any practical signal s(t) can be represented in the frequency domain by its Fourier 
transform S(f), given by 



A f°° 

S(f)=F{ S (t)}= s{t)e~^dt. (1-1.7) 

-** J —oo 

For convenience, the relation between s(t) and S(f) may be written “ s(t) < — -+ S(f) ” 

or simply “s(t) <-> The Fourier transform (FT) is in general complex; its 

magnitude is called the magnitude spectrum and its phase is called the phase 
spectrum. The square of the magnitude spectrum is the energy spectrum and 
shows how the energy of the signal is distributed over the frequency domain; the 
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total energy of the signal is 

/ oo /*oo 

\S(f)\ 2 df= S(f)S*(f)df (1.1.8) 

-oo J — oo 

where the superscripted asterisk (*) denotes the complex conjugate. Although the 
representation S(f) is a function of frequency only — time having been “integrated 
out” — the FT is a complete representation of the signal because the signal can be 
recovered from the FT by taking the inverse Fourier transform (IFT): 

/ OO 

S(f)e^df. (1.1.9) 

-oo 

But the “completeness” of the FT representation does not make it convenient for all 
purposes, as may be seen by considering the same three examples: 

1. Sinusoidal FM signal: If f> = ip = 0 in Eq. (1.1.4), the expression for s\(t) 
can be expanded into an infinite series as 

OO 

Slit) = A ^2 Jn{(3) cos27t(/ c +nf m )t (1.1.10) 

n =— oo 

where (3 is the modulation index (/3 = /d// m ) and J n denotes the Bessel func- 
tion of the first kind, of order n [6, p.226]. In the frequency domain, this 
becomes an infinite series of delta functions; one of these (the carrier) is at 
the mean frequency f c , and the remainder (the sideband tones) are separated 
from the carrier by multiples of the modulating frequency / m . Although the 
number of sideband tones is theoretically infinite, the significant ones 1 may 
be assumed to lie between the frequencies f c ± (/d + / m ) or, more conser- 
vatively, f c ± (/d + 2/ m ) . This information is essential if one is designing a 
tuning filter to isolate the TV audio carrier or separate one FM channel from 
adjacent channels. But it is inadequate if one is designing a modulator or 
demodulator, because its connection with the FM law is even more obscure 
than that of Eq. (1.1.4). 

2. The magnitude spectrum of the linear FM signal: Eq. (1.1.5) is shown 
on the right side of Fig. 1.1.2(a), for ip = 0, T = 64s, /o = 0.1Hz and 
a = (3/640) Hz/s. The graph shows that magnitude is significant in the band 
corresponding to the frequency sweep range (0.1 Hz < / < 0.4 Hz), and fur- 
thermore that the energy is mostly confined to that band. However, the graph 
fails to show that the frequency is increasing with time. In other words, the 

1 For more detail, see Carlson [6], pp. 220— 37. The above description considers only positive 

frequencies; similar comments apply to the negative frequencies. Theoretically, the lower sideband 
tones belonging to the positive-frequency carrier extend into the negative frequencies, while the 
corresponding sideband tones belonging to the negative-frequency carrier extend into the positive 
frequencies; but such “aliased” components are negligible if f c and are appropriately chosen. 
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magnitude spectrum tells us what frequencies are present in the signal, but 
not the “times of arrival” of those frequencies (the latter information, which 
we call time delay, is encoded in the phase spectrum). 

The shortcomings of the magnitude spectrum may be seen even more clearly 
in Fig. 1.1.2(b), which shows a signal whose magnitude spectrum (right) is 
identical to that of the linear FM signal in Fig. 1.1.2(a), but whose appearance 
in the time domain (left) is very different from a linear FM. 2 

3. Similarly, the musical performance has a magnitude spectrum which tells 
us what frequencies are present, but not when they are present; the latter 
information is again encoded in the phase. The magnitude spectrum may 
exhibit as many as 120 peaks corresponding to the notes of the chromatic 
scale in the audible frequency range, and the relative heights of those peaks 
may tell us something about the tonality of the music (or whether it is tonal 
at all), but the timing of the notes will not be represented in the magnitude 
spectrum and will not be obvious from the phase spectrum. 

These three examples show that the frequency-domain representation “hides” the 
information about timing, as S(f) does not mention the variable t. 

1.1.4 Joint Time-Frequency Representation 

As the conventional representations in the time domain or frequency domain are 
inadequate in the situations described above, an obvious solution is to seek a rep- 
resentation of the signal as a two-variable function or distribution whose domain 
is the two-dimensional (t, f) space. Its constant-f cross-section should show the 
frequency or frequencies present at time t, and its constant-/ cross-section should 
show the time or times at which frequency / is present. Such a representation is 
called a time-frequency representation (TFR) or time-frequency distribu- 
tion (TFD). 

As an illustration of what is desired, Fig. 1.1.3 shows one particular TFD of the 
linear FM signal of Eq. (1.1.5), for if) = 0, T = 64s, /o = 0.1 and a = (3/640) Hzs -1 . 
The TFD not only shows the start and stop times and the frequency range, but also 
clearly shows the variation in frequency with time. This variation may be described 
by a function fi(t), called the instantaneous frequency (IF). A signal may have 
more than one IF; for example, Fig. 1.1.4 shows a TFD of a sum of two linear 
FM signals, each of which has its own IF. 3 These IF features are not apparent in 
conventional signal representations. 

Non-stationary signals for which a TFD representation may be useful occur not 
only in broadcasting, seismic exploration and audio, from which our three exam- 

2 The signal in part (b) of Fig. 1.1.2 was obtained from that in part (a) by taking the FFT, 
setting the phase to zero, taking the inverse FFT and shifting the result in time. It is not the 
product of a sine function and a cosine function. 

3 N.B.: In Fig. 1.1.4 and in all subsequent graphs of TFDs, the labels on axes are similar to 
those in Fig. 1.1.3, but may be more difficult to read because of space constraints. 
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Fs=1 Hz N=65 
Time-res=1 






Fig. 1.1.3: A time- frequency representation of a linear FM signal [Eq. (1.1.5)] with duration 65 sam- 
ples, starting frequency 0.1 and finishing frequency 0.4 (sampling rate 1Hz). The time-domain repre- 
sentation appears on the left, and the magnitude spectrum at the bottom; this pattern is followed in 
all TFD graphs in Part I of this book. 




Fig. 1.1.4: A time- frequency representation of two linear FM signals with close parallel FM laws, 
duration 512 samples, frequency separation 0.04 (sampling rate 1Hz). 



pies are taken, but also in numerous other engineering and interdisciplinary fields 
such as telecommunications, radar, sonar, vibration analysis, speech processing and 
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medical diagnosis. Time-frequency signal processing (TFSP) is the processing 
of signals by means of TFDs. 

1.1.5 Desirable Characteristics of a TFD 

The use of a TFD for a particular purpose is inevitably based on particular as- 
sumptions concerning the properties of TFDs. To ensure that the conclusions of 
the analysis are sound, these assumptions must be identified and verified. Gen- 
eralizing somewhat, we may say that the desirable properties of a TFD are those 
on which the most likely applications of TFDs depend. The following applications 
illustrate some likely uses of TFDs: 

• Analyze the raw signal in the (t. f) domain so as to identify its characteristics 
such as time variation, frequency variation, number of components, relatives 
amplitudes, etc. 

• Separate the components from each other and from the background noise by 
filtering in the (t, /) domain. 

• Synthesize the filtered TFD in the time domain. 

• Analyze specific components separately: 

— Track the instantaneous amplitude; 

— Track the instantaneous frequency; 

— Track the instantaneous bandwidth (spread of energy about the IF). 

• Choose a mathematical model of the signal, showing clearly the significant 
characteristics, such as the IF. 

These applications can be carried out using a TFD with the following properties: 
1(a) The TFD is real (because energy is real). 

1(b) The integral of the TFD over the entire ( t , /) plane is the total energy of the 
signal (so that energy is conserved). 

1(c) The integral over a rectangular region of the (t, f) plane, corresponding to a 
finite bandwidth and finite time interval, is approximately the energy of the 
signal in that bandwidth over that interval, provided that the bandwidth and 
interval are sufficiently large (see Section 1.2.5 for details). 

2 For a nronocomponent FM signal, the peaks of the constant-time cross-sections 
of the TFD should give the IF law which describes the signal FM law. 

3 For a multicomponent FM signal, the dominant peaks of the TFD should 
reflect the components’ respective FM laws; and the TFD should resolve any 
close components, as in Fig. 1.1.4. 
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A simple measure of the resolution of a TFD is the concentration of energy about 
the IF law(s) of a signal [Article 7.4 and Chapter 10]. 

The linear FM signal, because of its simple FM law, is a convenient test signal 
for verifying the first four properties. Property 3 may be tested by a sum of linear 
FM signals. 

The above properties have further implications. For example, the resolution 
property (3) helps to ensure robustness in the presence of noise, while the last two 
properties (2 and 3) together allow discernment of the multiple IFs of multicompo- 
nent signals. 

1.2 Signal Formulations and Characteristics in the (t, f) Domain 
1.2.1 Signal Models used in (t,f) Methods 

To represent signals such as the linear FM, several types of signal models are com- 
monly used in TFSAP. The choice of model depends on the number and nature 
of the parameters needed to describe the signal. For example, a single sinusoid 
with constant frequency and normalized amplitude and phase is described by the 
equation 

S 4 (t) = cos 27 rf c t , ( 1 . 2 . 1 ) 

in which the only parameter is the frequency f c . If the amplitude and phase are 
significant in the application, two more parameters are needed. A linear combination 
of such signals may be written in the form 

M 

Ss(t) = ^2 a k cos(27T fkt + ipk) , (1.2.2) 

k=l 



which involves 3 M parameters. 

Because S 4 (t) and 55 (f) contain terms (or “components”) of constant amplitude, 
frequency and phase, they are completely and clearly described by Fourier trans- 
forms; no time dependence is required. But a sinusoidal FM signal or chirp signal 
requires a TFD. A musical audio signal also calls for a form of TFD (as is suggested 
by the notation in which music is written) , and the TFD should clearly resolve the 
multiple components. 

Further difficulties are raised by more complex signals such as 



se(t) 



M 



^2a k (t)eP 2 ^ o/*(r)d r + w (t) 



\k = 1 



(1.2.3) 



where Ofe(f) is the time-varying amplitude of the k th component, fk(t) is the time- 
varying frequency of the k th component, and w(t) is additive noise. The analysis of 
such a signal requires the ability not only to distinguish the time-varying compo- 
nents from each other in spite of their varying amplitudes and frequencies, but also 




Signal Formulations and Characteristics in the ( 7 , /) Domain 



13 



to separate them from the noise. Such comments still apply if the amplitude of the 
fc th component is a multiplicative noise factor nik{t), as in the signal 

/ m \ 



s r {t) = mk (*) e32 " J ° fk{T)dT + (1.2.4) 



\fe=i / 

In this case, an even more specialized analysis will be required. 



1.2.2 Analytic Signals 

It is well known that a signal s(t) is real if and only if (iff) 

= (1.2.5) 

where S(f) is the Fourier transform of s(t). In other words, a real signal is one 
that exhibits Hermitian symmetry between the positive-frequency and negative- 
frequency components, allowing the latter to be deduced from the former. Hence the 
negative-frequency components of a real signal may be eliminated from the signal 
representation without losing information. In the case of a real lowpass signal, 
removal of negative frequencies has two beneficial effects. First, it halves the total 
bandwidth, allowing the signal to be sampled at half the usual Nyquist rate without 
aliasing [7,8]. Second, it avoids the appearance of some interference terms generated 
by the interaction of positive and negative components in quadratic TFDs (to be 
treated in detail in Section 3.1.2). 

Definition 1.2.1: A signal z(t ) is said to be analytic iff 

Z(f) =0 for f< 0 , (1.2.6) 

where Z(f) is the Fourier transform of z(t). 

In other words, an analytic signal contains no negative frequencies', it may have a 
spectral component at zero frequency (DC). 

Theorem 1.2.1: The signal 

z(t) = s(t) + jy{t) , (1.2.7) 

where s(t) and y(t) are real, is analytic with a real DC component, if and only if 

Y(f) = (-jsgnf)S(f) (1.2.8) 

where S(f) and Y(f) are the FTs of s(t) and y(t), respectively, and where 

. f -1 if Z < 0 ; 

sgn^ = ^ 0 if (, = 0; (1.2.9) 

l+i */ £ > o . 



Proof: Take the FT of Eq. (1.2.7) and use Eq. (1.2.6). 
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1.2.3 Hilbert Transform; Analytic Associate 

If the Fourier transforms of s(t) and y(t) are related according to Eq. (1.2.8), we 
say that y(t) is the Hilbert transform of s(t), and we write 

y(t) = H{s(t)}. (1.2.10) 

Hence we may restate Theorem 1.2.1 as follows: A signal is analytic with a real DC 
component if and only if its imaginary part is the Hilbei't transform of its real part. 

By invoking the “if” form of Theorem 1.2.1 and restating the sufficient condition 
in terms of Eq. (1.2.7), we may now see the practical significance of the theorem 
and the Hilbert transform: Given a real signal s(t), we can construct the complex 
signal 

z(t) = s(t) + jH{s(t)} (1.2.11) 

and know that z(t) is analytic. This z(t) is called the analytic signal “corresponding 
to” or “associated with” the real signal s(t). In this book, for convenience, we shall 
usually call z(t) the analytic associate of s(t). 

By taking the IFT of Eq. (1.2.8) and applying Eq. (1.2.10), we arrive at the 
following concise definition of the Hilbert transform: 

Definition 1.2.2: The Hilbert transform of a signal s(t), denoted by TL{s(t)}, 

is 

H{s(t)}=T {(-jsgn /)£{«(*)}}. (1.2.12) 

where T {. . .} denotes the Fourier transform. 

In other words, the Hilbert transform of s(t) is evaluated as follows: 

1. Take the Fourier transform S(f) of s(f); 

2. Multiply S(f) by —j for positive /, by +j for negative /, and by zero for 

/ = 0; 

3. Take the inverse Fourier transform. 

According to step 2 of the above procedure, a Hilbert transformer introduces a 
phase lag of 90 degrees (as — j = e - -' 71 '/ 2 ), producing a signal in quadrature to the 
input signal. The effect is well illustrated by the following result, which is easily 
verified using Definition 1.2.2 and a table of transforms: 

Example 1.2.1: If fo is a positive constant, then 



7t{cos(27r/ 0 t)} = sin(27r/ 0 t) 
7d{sin(27r/ 0 t)} = - cos( 27 t f 0 t) 



(1.2.13) 

(1.2.14) 
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It would be convenient if the pattern of Example 1.2.1 were applicable to mod- 
ulated signals; so that for example, we could say 

H{a(t) cos</>(f)} = a(t) sin<fi(t) , (1.2.15) 

which would imply that the analytic associate of the real signal s(t) = a(t) cos <j>(t) 
is 



z(t) = a(t) cos (j>(t) + jTi{a(t) cos cj>(t)} 

= a(t) cos cj)(t) + ja(t) sin <j>(t) 

= a(t)e m) . (1.2.16) 

The condition under which Eq. (1.2.15) holds is that the variation of a(t) is suf- 
ficiently slow to ensure “spectral disjointness”, i.e. to avoid overlap between the 
spectrum of a(t) and the spectrum of cos 

Eq. (1.2.8) indicates that the transfer function of a Hilbert transformer is 
— jsgaf. The corresponding impulse response is 

r~U-jsgnf} = - (1.2.17) 

*<-/ 7 rt 



Using this result and applying the convolution property to Eq. (1.2.12), we obtain 
a definition of the Hilbert transform in the time domain: 



H{s(t)} 



s(t) , - 




(1.2.18) 

(1.2.19) 



where p.v.{. . .} denotes the Cauchy principal value of the improper integral [9], 
given in this case by 4 



lim 

<5-*0 




s(r) 

t — T 



dr + 




( 1 . 2 . 20 ) 



1.2.4 Duration, Bandwidth, BT Product 
1.2. 4.1 Finite Duration and Finite Bandwidth 

In real life, signals are observed for finite periods of time and processed by devices 
with finite usable bandwidths. If a signal s(t) has the Fourier transform £(/), the 
duration of the signal is the smallest range of times outside which s(t) = 0 , while 
the bandwidth of the signal is the smallest range of frequencies outside which 
S(f) = 0. These definitions, as we shall see, lead to the conclusion that a finite 

4 In practice Eq. (1.2.18) is rarely used, because the FT properties make it easier to work in the 
frequency domain. 
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duration implies infinite bandwidth and vice versa. In practice, however, any useful 
signal must have a beginning and an end (finite duration) and its FT must be within 
the frequency range of the measuring/processing equipment (finite bandwidth). 5 In 
practice, therefore, the strict definitions need to be relaxed in some way. 

A time-limited signal, of duration T centered at time t = 0 , can be expressed 
as 

s T (t) = s(t) rect[f/T] , (1.2.21) 

where the subscript T indicates the duration [see Eq. (1.1.6)]. 

The FT of s T (t) is 

S T (f) = S(f) * T sine fT (1.2.22) 

where * denotes convolution in frequency. Thus the bandwidth of St { /) is infinite. 

If, in order to avoid the effects of discontinuities, we replace rect[f/T] with a 
smoother window w(t) of the same duration T, we can write 

s T (t) = s(t)rect[t/T]w(t) , (1.2.23) 

whose Fourier transform still involves a convolution with sine fT, giving an infinite 
bandwidth. 

In analogy to the time-limited case, a band-limited signal, of bandwidth B 
centered at the origin, can be expressed in the frequency domain as 

S B (f) = S(f)rect{f/B\ (1.2.24) 

In the time domain, the signal is given by 

s B (t) = s(t) * B sine Bt (1.2.25) 

which has an infinite duration. Thus, under the “usual” definitions of duration and 
bandwidth, a finite bandwidth implies infinite duration. 

1.2. 4. 2 Effective Bandwidth and Effective Duration 

If there is no finite bandwidth containing all the energy of the signal, there may 
still be a finite bandwidth containing most of the energy. Hence, for example, 
a bandwidth containing 99% of the signal energy might be accepted as a useful 
measure of the signal bandwidth. If the nominated fraction of the signal energy 
were confined between the frequencies / m i n and / max , the bandwidth would be 

B — ,/max fmin • 

A less arbitrary but more relaxed measure of bandwidth is the so-called effective 
bandwidth B e , defined by 

i r°° 

B *=Tf\ /W)| 2 # (1-2.26) 

J — oo 



5 In practice, all acquisition and measuring systems are low-pass. 
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where S(f) is the FT of the signal, and E s is the total energy of the signal, given 
by 



E s = \S(f)\ 2 df. 



(1.2.27) 



B 2 is the second moment of |£(/)| 2 w.r.t. frequency, about the origin (/ = 0). 
For brevity, we call B 2 the “second moment of the signal w.r.t. frequency”. As an 
aid to remembering the definition, note that if / were a random variable and |<!5(/)| 2 
were its probability density function (p.cl.f.), then we would have E s = 1, so that B 2 
would be the variance of / if the mean of / were zero. Thus the effective bandwidth 
B e is analogous to the standard deviation of /. When we consider how little of a 
typical probability distribution falls within one standard deviation, we realize that 
the effective bandwidth is only a mathematical construction and should not be used 
as an estimate of the bandwidth required for accurate measurement and processing 
of the signal. 

The definition of duration, like that of bandwidth, can be relaxed so as to ob- 
tain both finite bandwidth and finite duration. For example, one could define the 
duration T where a nominated fraction of the signal energy occurs. The so-called 

effective duration T e , defined by 



Tg 



1 

El 




t 2 \s(t)\ 2 dt , 



(1.2.28) 



is the second moment of |s(t)| 2 w.r.t. time, about the origin (t = 0). For brevity, 
we refer to T 2 as “the second moment of the signal w.r.t. time.” The definitions of 
effective duration and effective bandwidth were proposed in 1946 by Gabor [10]. 6 

Slepian [11] has proposed definitions of duration and bandwidth based on the ac- 
curacy of the detecting and measuring apparatus. Let us define a “time truncation” 
of s(t) as a signal s(t) that satisfies 



[ 0 if t < fi 

s(t) = < s(t) if t,\ < t < <2 (1.2.29) 

( 0 if t > t-2 

where G > t \ , so that the duration of s{t) is t? — . Let a “frequency truncation” be 

defined similarly in the frequency domain. Then the Slepian duration of s(f) is the 
duration of the shortest time-truncation s(t) that the apparatus cannot distinguish 
from s(t), while the Slepian bandwidth of s(t) is the bandwidth of the most 
narrow-band frequency-truncation S(f) that the apparatus cannot distinguish from 

S(f). 

Given suitable definitions of B and T, the bandwidth-duration product BT 
is self-explanatory. The BT product also has a practical significance: because a 

6 Dennis Gabor (1900-1979), a Hungarian-born electrical engineer who settled in Britain, is 
best known for the invention of holography (1947-8), for which he was awarded the Nobel Prize 
for Physics in 1971. 
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signal of total bandwidth B can be reconstructed from samples at the sampling rate 
B, the total number of samples required to represent the signal — i.e. the number 
of degrees of freedom in the signal— is equal to BT. In other words, BT is a 
measure of the information richness of the signal. This implies that there may be 
little point in attempting to extract large amounts of information from a signal with 
a small BT product. 



Example 1.2.2: For the Gaussian signal 



s(t)= e - a2t \ 

the effective duration is 

T e = — . 
2a 


(1.2.30) 


(1.2.31) 


The FT of s(t) is 




s(f) = ^ e -* 2 f 2 / a \ 

a 


(1.2.32) 


so that the effective bandwidth is 




B e = — . 
2tt 


(1.2.33) 


From Eqs. (1.2.31) and (1.2.33), 




B e T e = - — • 
4? r 


(1.2.34) 



It can be shown that the Gaussian signal is the only signal for which this equality 
holds, and that for all other signals B e T e > [10]. 



1.2.5 Asymptotic Signals 

For signals with significant information content, it is desired to know not only the 
overall bandwidth B , but also the distribution of energy through the bandwidth, e.g. 
the frequencies present, their relative amplitudes, and the times during which they 
are significant. Similarly we may want to know not only the overall duration T, but 
also the distribution of energy throughout the duration, e.g. the times during which 
the signal is present, the relative amplitudes at those times, and the significant 
frequencies present during those times. Such signals may be modeled by the class 
of asymptotic signals. 

Definition 1.2.3: A signal s(t) is asymptotic iff it has the following properties: 

(a) The duration T , as defined by Slepian, is finite; 

(b) The bandwidth B, as defined by Slepian, is finite; 

(c) The product BT is large (e.g. greater than 10); 
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(d) The amplitude is bounded so that 




is finite. 



(1.2.35) 



Asymptotic signals allow useful approximations for deriving simple signal models 
(e.g. to express analytic signals). 



1.2.6 Monocomponent vs. Multicomponent Signals 

A monocomponent signal is described in the (t, /) domain by one single “ridge”, 
corresponding to an elongated region of energy concentration. Furthermore, inter- 
preting the crest of the “ridge” as a graph of IF vs. time, we require the IF of a 
nronocomponent signal to be a single-valued function of time. Fig. 1.1.3 shows an 
example of a nronocomponent signal. 

Such a nronocomponent signal has an analytic associate of the form 

z(t) = aft) eP^' 1 , (1.2.36) 

where (fit) is differentiable, and aft) is real and positive (being the instantaneous 
amplitude). If sft) itself is real and asymptotic, it can be expressed as 

s(t) = aft) cos (fft). (1.2.37) 

A multicomponent signal may be described as the sum of two or more nronocom- 
ponent signals. Fig. 1.1.4 shows an example of a multicomponent signal (composed 
of two components). 



1.3 Instantaneous Frequency and Time-Delay 
1.3.1 Instantaneous Frequency (IF) 

Definition 1.3.1: The instantaneous frequency of a monocomponent signal is 

m = ^c/>'(t) (i.3.i) 

where (fft) is the instantaneous phase of the signal. 

This formulation is justified below by considering first a constant-frequency signal, 
then a variable-frequency signal. 

Consider the amplitude-modulated signal 

xft) = aft) cos ( 27 t f c t + if) ( 1 . 3 . 2 ) 

where f c and if are constant. As t increases by the increment 1 / f c , the argument of 
the cosine function increases by 27T and the signal passes through one cycle. So the 
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period of the signal is l// c , and the frequency, being the reciprocal of the period, 
is f c . The same signal can be written as 



where 



from which we obtain 



xft) = aft) cos (f>(t) 


(1.3.3) 


(f{t) = 2ir f c t + if , 


(1.3.4) 




(1.3.5) 



Although the left-hand side of this equation (the frequency) has been assumed 
constant, the right-hand side would be variable if <j>(t) were a nonlinear function. 
So let us check whether this result can be extended to a variable frequency. 
Consider a signal whose analytic associate is 



zft) = aft) 



(1.3.6) 



where a(t) and <j>(t) are real and a(t) is positive; then aft) is called the instan- 
taneous amplitude and <fi(t) is called the instantaneous phase. Let <j>(t ) be 
evaluated at t = t\ and t = t 2 , where t -2 > t\ . By the mean value theorem of el- 
ementary calculus, if (fft) is differentiable, there exists a time instant t between t\ 
and t '2 such that 

4>(t 2 ) - (ffti) = ft 2 - ii) (j>'ft). (1.3.7) 

Let pi be the period of one particular oscillation of z(t), and let /) = 1 /p, . If 
t 2 =ti+Pi, then (j>ft 2 ) = <j>(ti) + 27 r, so that Eq. (1.3.7) becomes 



27T = Pif> ft) ; 



(1.3.8) 



that is, 



, _ M 

h 2t r 



(1.3.9) 



Now t is an instant during a cycle of oscillation and /) is the frequency of that 
oscillation, suggesting that the right-hand side be defined as the instantaneous 
frequency (IF) at time t, as in Definition 1.3.1 above. 

Comparing Eqs. (1.2.37) and (1.3.6), we see that 



s(t) = Re{ 2 (t)} . 

where Re{. . .} denotes the real part. Now let us define 

yft) = Im-{+(f)} = a(t) sin c/)(t) 



(1.3.10) 



(1.3.11) 



where Im{. . .} denotes the imaginary part. 

Using Eq. (1.3.1), we can easily confirm that the signals described by Eqs. (1.1.4) 
and (1.1.5) have the IFs given in Eqs. (1.1.1) and (1.1.2) respectively. 
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Definition 1.3.1 is strictly meaningful only for a monocomponent signal, while 
a multicomponent signal ought to have a separate IF for each component. In par- 
ticular, note that the IF of the sum of two signals is not the sum of their two 
IFs. 

The IF is a detailed description of the frequency characteristics of a signal. This 
contrasts with the notion of “average frequency” defined next. 



Definition 1.3.2: The average frequency of a signal is 



fo 



f-OO 

/ f\S(f)\ 2 df 
Ju 

/ \S(f)\ 2 df 

Jo 



where S(f) is the Fourier transform of the signal. 



(1.3.12) 



In other words, /o is the first moment of \S(f)\ 2 w.r.t. frequency; that is, we 
define the “average” frequency as if | S(f) \ 2 were the probability density function of 
the frequency. 

Notice that if | S(f)\ 2 were replaced by a TFD, f 0 would become a function 
of time, suggesting that perhaps the first moment of a TFD w.r.t. frequency is 
a measure of IF. The conditions under which this is true will be explained in due 
course (and summarized in the table on p. 75). For the moment, we simply note that 
any reasonable time- frequency representation of a signal should contain information 
about the IF laws of the components. In particular, it would be most convenient if 
the crests of the ridges in the (t, /) domain represented the IF laws. 



1.3.2 IF and Time Delay (TD) 

The IF of a signal indicates the dominant frequency of the signal at a given time. 
Let us now seek a dual or “inverse” of the IF, indicating the dominant time when 
a given frequency occurs. 

If z{t) is an analytic signal with the Fourier transform 

Z(f) = A5(f — fi) (1.3.13) 

where A is in general complex, the dominant frequency is fi. 

Taking the inverse FT of Z(f) gives 

z[t) = Ae j2nfit . (1.3.14) 

The instantaneous phase of z(t), denoted by <j>(t), is 

( t>(t ) = arg z(t) = fit + arg A (1.3.15) 



fi 




so that 



(1.3.16) 
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Although this result has been obtained for constant frequency, the right-hand of 
Eq. (1.3.16) is also valid for a variable-frequency signal, as explained earlier. 

Now let us repeat the argument with the time and frequency variables inter- 
changed. The signal 

z(t) = a6(t — Td) , (1.3.17) 

where a is in general complex, is an impulse at time Td. If we ask what is the “delay” 
of this signal, the only sensible answer is Td . The FT of this signal is 

Z(t) = ae- j2 * fTA (1.3.18) 

and its phase, denoted by 0(f), is 

0(f) = arg Z(f) = -27r/T d +arga (1.3.19) 

so that 

Td = ~h 9 ' if) - (L3 ' 20) 

Again, although this result has been obtained for a constant Td , the right-hand side 
of Eq. (1.3.20) is well defined even if 6'(f) varies with f. 

Definition 1.3.3: If z(t) is an analytic signal with the Fourier transform Z(f), 
then the time delay (TD) of z(t), denoted by Td (/), is 

Td(f) = ~^0'(f) , (1.3.21) 



where 0(f) = arg Z(f) 

Notice that the definitions of Td and ,/) are similar, except that time and fre- 
quency are interchanged and Eq. (1.3.21) has an extra minus sign; hence we say that 
time delay is the dual of instantaneous frequency. 7 

Seeing that the instantaneous frequency fi(t) is a function assigning a frequency 
to a given time, whereas the time delay Td(/) is a function assigning a time to a 
given frequency, we may well ask whether the two functions are inverses of each 
other. Clearly they are not always inverses, because the IF function may not be 
invertible. So let us restrict the question to invertible signals, i.e. monocomponent 
signals whose IFs are monotonic functions of time. 

One example of an invertible signal is the generalized Gaussian signal, i.e. a linear 
FM signal with a Gaussian envelope. Let such a signal peak at time to, with peak 
amplitude A, center frequency f c , sweep rate a and decay constant (3, and suppose 

7 The term “time delay” is well established, but tautological. The term “instantaneous frequency” 
is a quantity with the dimensions of frequency modified by an adjective indicating localization in 
time. The dual of this term should therefore be a quantity with the dimensions of time modi- 
fied by an adjective indicating localization in frequency, e.g. “frequency-localized delay” or, more 
succinctly, “frequency delay”. 
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that the instantaneous frequency is positive while the envelope is significant. Then 
the analytic associate of the signal is 

z(t) = A exp (j27r f c [t - t 0 ] + ^ 4 ^% - *o] 2 ) (1.3.22) 

and using Eq. (1.3.16) with a(t) = A exp (— n/3(t — to) 2 ), its IF is 

Mt) = f c + a[t - t 0 \. (1.3.23) 

To find the inverse function of fi(t), we simply solve for t, obtaining 

t = t 0 + ^ ~ , (1.3.24) 

a 

which suggests that the time delay 7d(/) of z(t) can be estimated by 

Td(f) = to + (1.3.25) 

a 

In general 7d and ra are not equal, but converge as the BT product increases. As an 
example, Fig. 1.3.1 shows the IF (solid line) and TD (dotted line) of the generalized 
Gaussian signal z(t) for two values of B e T e , where the subscript e means “effective”. 
Note that the two curves are closer for the larger B e T e . 

1.3.3 Mean IF and Group Delay (GD) 

Let z(t) be a bandpass analytic signal with center frequency / c . Let its Fourier 
transform be 

Z(f) = M(f-f c )e^ (1.3.26) 

where the magnitude M(f — f c ) and phase 6(f) are real. If the signal has linear 
phase in the support of Z(f), i.e. if 6(f) is a linear function of / wherever Z(f) is 
non-zero, we can let 

0(f) = — 27TTp/ c - 27TTg[/ - / c ] (1.3.27) 

where r p and r g are real constants with the dimensions of time. Eq. (1.3.26) then 
becomes 

Z(f) = M(f — f c ) e -J'( 27rT p/ c + 27rT g[/-/c]) (1.3.28) 

= e ~ j2 ^ f ^ M(f - f c ) e - j2nT ^ f - f J. (1.3.29) 

Taking the IFT of Z(f ), we find 

z(t) = m(t - r g ) e J ' 27r A[t-T P ] (1.3.30) 

where m(t) is the IFT of M(f). Now because M(f) is real, m(t) is Hermitian 
[i.e. m(—t) = m*(t )], so that |m(f)| is even. Hence r g is the time about which 
the envelope function is symmetrical; for this reason, r g is called the group de- 
lay [12, pp. 123-124]. The phase of the oscillatory factor is — 2irf c T p , wherefore r p 
is called the phase delay. These observations together with Eq. (1.3.27) lead to 
the following definition: 
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IF (solid) and Group Delay (dotted) 




IF (solid) and Group Delay (dotted) 




Fig. 1.3.1: Instantaneous frequency and time delay for a linear FM signal with a Gaussian envelope: 
(a) total duration T — 33 secs, B e T e = 0.1806; (b) T = 65 secs, B e T e = 0.3338. For each signal, 
the left-hand graph shows the time trace, spectrum and TFD [Wigner-Ville distribution, defined later 
in Eq. (2.1.17)], while the right-hand graph shows the IF (solid line) and time delay (dotted line). The 
vertical dotted segments are caused by truncation of the frequency range (to avoid finite-precision 
effects). 



Definition 1.3.4: If an analytic signal has the Fourier transform, 



Z(f)=\Z(f)\e je W , 

then the group delay (GD) of the signal is 

T g (/)=~^'(/) 



( 1 . 3 . 31 ) 



( 1 . 3 . 32 ) 
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and the phase delay of the signal is 

T P (f) = ^^jO(f). (1.3.33) 

Eq. (1.3.32) is found by differentiating Eq. (1.3.27) w.r.t. /, and Eq. (1.3.33) is found 
by putting f c = f in Eq. (1.3.27). Whereas Eq. (1.3.27) assumes linear phase, the 
above definition is meaningful whether the phase is linear or not. 

Further, Eq. (1.3.32) is the same as Eq. (1.3.21). But the physical interpretations 
are different; the time delay applies to an impulse, whereas the group delay applies 
to the envelope of a narrowband signal. 

Now consider the dual of the above argument. Let zft) be a time-limited signal 
centered on t = t c , and let 

zft) = aft — t c ) e iW) (1.3.34) 

where aft) and cf>ft) are real. If the signal has constant instantaneous frequency 
in the support of zft), i.e. if <j>ft) is a linear function of t wherever zft) is non-zero, 
we can let 

<t>ft) = 27r/ 0 f c + 2 t Tf m [t - t c ] (1.3.35) 

where /o and f m are real constants with the dimensions of frequency. Eq. (1.3.34) 
then becomes 

zft) = a ft — t c ) e K 27r /o*c+27r/ m [t-t c ]) (1.3.36) 

= e^^aft-Qe^^^l (1.3.37) 

Taking the FT of zft), we find 

Zff) = Aff - f m ) e -^ 2 -[/-/o]* c (4.3.38) 

where Aff) is the FT of aft). Now because aft) is real, Aff) is Hermitian, so that 
|H(/)| is even. Hence f m is the frequency about which the amplitude spectrum is 
symmetrical; for this reason, f m is called the mean IF. Differentiating Eq. (1.3.35) 
w.r.t. t, leads to the following definition: 

Definition 1.3.5: For the signal 

zft) = \zft)\e m) , (1.3.39) 

the mean IF is 

fmft) = ^<f'ft). (1.3.40) 

Thus the mean IF is the same as the IF defined earlier [Eq. (1.3.1)], but the 
physical interpretations are different. The IF has been derived for a tone (and 
earlier for a modulated sinusoid), whereas the mean IF applies to the spectrum of 
a short-duration signal. 
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1.3.4 Relaxation Time, Dynamic Bandwidth 

For a linear FM signal, the instantaneous phase <j>{t) is quadratic. So <j>(t) can be 
expanded in a Taylor series about t = to- 

4>{t) = + </>'(t 0 )[t-t 0 ] + \(i>" {t 0 )[t-t 0 } 2 (1.3.41) 

= </)(to) + 2 TTfi(t 0 )[t-t 0 ] + ^2nfi(to)[t—to] 2 . (1.3.42) 

The relaxation time T r , as defined by Rihaczek [13, p. 374], is the duration over 
which the instantaneous phase deviates no more than 7t/4 from linearity. That is, 

|§27r/ i / (i 0 )T I ?/4|=7r/4. (1.3.43) 

Solving this equation leads to the following definition: 

Definition 1.3.6: The relaxation time of a signal is 

T r (t)= ' (1.3.44) 

at 

where fi(t ) is the instantaneous frequency. 

The dual of relaxation time, known as dynamic bandwidth, is the bandwidth 
over which the phase spectrum, assumed to be a quadratic function of frequency, 
deviates no more than 7t/4 from linearity. The result is as follows: 

Definition 1.3.7: The dynamic bandwidth of a signal is 

BM ) = ^ ^ (1.3.45) 

df 

where t<j(/) is the time delay. 

As the relaxation time is a measure of the time needed to observe significant 
variations in IF, so the dynamic bandwidth is a measure of the bandwidth needed 
to observe significant variations in time delay. 

1.4 Summary and Discussion 

Clear rationales are developed for justifying the need to use joint time-frequency 
representations for non-stationary signals such as FM signals. 

An understanding of the concept of instantaneous frequency and its dual, time 
delay, is necessary for the interpretation of TFDs. The use of the analytic associate 
of a given real signal, rather than the signal itself, is useful for reducing the required 
sampling rate, and essential for obtaining an unambiguous instantaneous frequency. 

The analytic associate is obtained via the Hilbert transform. The notions of BT 
product and monocomponent signals are introduced. For monocomponent asymp- 
totic signals, analytic associates can be written using an exponential form. The 
next chapter introduces various formulations of time-frequency distributions and 
demonstrates the importance of using the analytic associate. 
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